Glauber Critical Dynamics: Exact Solution of the Kinetic Gaussian Model 
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In this paper, we have exactly solved Glauber critical dynamics of the Gaussian model on three 
dimensions. Of course, it is much easy to apply to low dimensional case. The key steps are that we 
generalize the spin change mechanism from Glauber's single-spin flipping to single-spin transition 
and give a normalized version of the transition probability . We have also investigated the dynamical 
critical exponent and found surprisingly that the dynamical critical exponent is highly universal 
which refer to that for one- two- and three-dimensions they have same value independent of spatial 
dimensionality in contrast to static (equilibrium) critical exponents. 
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I. INTRODUCTION 



Irreversible dynamic systems exhibit complicated and interesting non-equilibrium phenomena near the critical point. 
The study of non-equilibrium statistical mechanics is much more difficult than equilibrium state due to the complexity. 
However, the interesting dynamic critical behaviors have been attracting a large number of researchers to work hard 
S3 ■ for many years. 

| Up to now, there is no general theory based on the first principle to describe the dynamic problems. However, a 
great progress has been made since the pioneering work completed by Glauber Jl]] and Kawasaki ||. According to 
their theory, the time-evolving of the order parameters is described by Markov processes with Glauber single spin 
O , flipping mechanism or Kawasaki exchange mechanism between two spins. Since then people have been paid much 
attentions to the study of critical dynamics. The research so far have been extended from the kinetic Ising model to 
t—{ | the kinetic Potts model and from integer to fraction dimension, in which many approximate methods such as Monte 
^ ■ Carlo simulation, high-temperature series expansion, e-expansion, bond-moving rcnormalization-group method, etc. 
\D ; have been applied 

Now let us turn on master equation, a basic equation for treating critical dynamics. As we know, the key step for 
solving master equation is the determination of the transition probability. Usually the transition probability between 
different states, i.e., different spin-configurations of the system, is only determined in terms of the detailed balance 

■ condition. Since such a choice is not unique, then leaves behind some arbitrariness. For removing the arbitrariness, 
at least in part, we suggest a normalized transition probability, this means in unit time interval the transition may 

■ occurs or may not. We apply this point to continuous spin (Gaussian) system and obtain exact solutions of one-, two- 
and three dimensional kinetic Gaussian model. 

This paper is organized as follows. In Sec.[n], a single-spin transiting critical dyn amics which suits arbitrary spin 
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systems is presented. As an application, the Gaussian model is treated in Sec. III. We not only obtain the exact 



solutions of ID, 2D and 3D kinetic Gaussian model, but also find that the dynamic critical exponent is highly 



universal. Finally, the conclusion is devoted to Sec. IV 



II. FORMALIZATION 



For an irreversible dynamic system subjected to the time-dependent perturbation, once the perturbation is removed, 
the system will very slowly approach to equilibrium state because of the large-scale fluctuation near the critical point, 
which is what we call the critical slowing down phenomenon. The critical dynamics will attempt to explain why the 
short-range interactions lead to long-time relaxation. Due to the complexity there have been no microscopic theory 
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based on the first principle so far, thus a suitable theoretical model will be quite important. As already mentioned, 
both the Glauber dynamics and Kawasaki dynamics have proven to be successful in many dynamic systems. 

In this section, we give a brief review of Glauber dynamics. Of course, we will give some improvement so that it can 
suits arbitrary-spin systems. For clearness, we start from the one-dimensional case (ID), and then the formulation is 
easy to extend to 2D and 3D. 

The ID lattice-spin model we shall discuss is a stochastic one. The spins of N fixed particles are represented as 
stochastic functions of time &j(t), (J — 1, N), which can be taken discrete values (discrete-spin model) or continuous 
values (continuous-spin model), and make transitions among these values. The transition, according to the Glauber dy- 
namics, can only change single spin value each time, such as <Tj (t) — > <Zj (t) because of the interacting of the system with 
heat reservoir. The transition probability Wj(<Tj(t) — > Oj(t)) from configuration (ai(t), 02(i), • • • , Cj(t), ■ ■ ■ , crjv(t)) to 
configuration (ui(<), oaCO, ' ' ' > ^.j W> ' ' ' > a N(t)), in general, depends on the momentary values of the neighboring spins 
as well as on the influence of the heat bath. For this reason statistical correlations exist between different spins. There- 
fore, it is necessary to deal with the entire TV-spin system as a unit. The evolution of spin functions describing system 
form a Markov process of N discrete or continuous random variables with a continuous time variable as argument. 

We introduce a probability distribution function p(cri, • • • , erjv, t), which denotes the probability of spin system being 
in the state (cxi, • • • ,<x/v) at time t. Let Wj((Tj — > &j) be the transition probability per unit time that the j'th spin 
transits from one value <Xj to another possible value &j, while the others remain fixed. Then, on the supposition of 
single-spin transition, we may write the time derivative of the function P{cri, • • ■ , (?n , t) as 

^P({o},t) XX 11 a MMit) + Wji&j - <7 j )p({<r w },ff il i)] . (I) 

3 &j 

This is a probability equation, in which the first term in the hand-right side of Eq.([j]) denotes the decrease of the 
probability distribution function P({a},t) per unit time due to the transition of the spin state from the initial value 
Uj (J — 1, 2, • • • , N) to various possible final values &j, and the second term denotes the increase of the probability 
distribution function P({a}, t) per unit time due to the transition of the spin state from the various possible initial 
values frj, (J = 1, 2, • • • , N) to final value <jj. We shall refer to the Eq.(|]) as the master equation since its solution 
would contain the most complete description of the system available. 

It is the most crucial step, obviously, that the transition probability must be determined before the master equation 
can be solved. Then, how to determine the transition probability? For this problem, Glauber' theory leaves some 
leeway to select. However, inappropriate selection will probably make the problem difficult to solve. So, we hope 
to find out a more definite expression to apply the Glauber' theory to arbitrary-spin systems. Now we consider it 
from both mathematics and physics aspects. In mathematics, generally speaking, the probability must be ergodic 
and positive definite, and can be normalized; in physics, we often request that an equilibrium thermodynamic system 
satisfies the detailed balance condition. Based on these considerations, we can restrict the spin transition probability 
Wj (<jj — > )tA to satisfy the following conditions so that ensure the system in thermodynamic equilibrium state 
finally: 

For the set (S, S) composed of a subset S and its dual subset S in phase space, existing <7j belonging to S and &j 
belonging to S , we have 



(1) ergodicity 

(2) positivity 

(3) normalization 

(4) detailed balance 



Va 3 ,a r . Wfci ->*,•) ?4 0; (2) 



Vo-,-,3*: W^fo-S^O; (3) 



Va, : £Vi fa -►*,■) = 1; (4) 



Vcr . & . . W^CTj -> Oj) = P eq (cri, ■■■;%) ,^ 
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in which 



P eq = | exp[-/3W(M)], Z = Y, exp[-/3H({a})], 

M 

where P eg is the equilibrium Boltzmann distribution function, Z is the partition function, 7i({er}) is the system 
Hamiltonian. 

Although the spin transition probabilities are not determined uniquely by above restriction conditions, there are 
less room to choose them. Furthermore, considering the fact that the transition of the individual spin depends merely 
on the momentary values of the neighboring spins as well as on the influence of the heat bath, we can suppose that 
the transition probability from <Tj to &j depend only on the heat Boltzmann factor of the neighboring spins, i.e. 



Wiipi — > &C) oc exp 



-(3Hi t,. a 3 



<*j> 



or 



Wi((Ji -> &i) = — exp 



-/3Hi [ &i, ^2 a i 
<ij> 



(6) 



where X)<ij> means that the summation for j is only related to the neighboring values of i. By means of the 
normalized condition (Eh, the normalized factor Qi can be determined as 



exp 



(7) 



Obviously, is independent of <7j and is related to the temperature and neighboring spins. 

Compared with Glauber' expression[l], Eq.(^) is a normalized version of transition probability. As we know, the 
constant a in Glauber' expression is a free constant determined by the time scale. Actuaclly, our expression is only a 
definite selection for constant a by extra restriction conditions and physical considerations. 

Usually, we are interested in local magnetization and the spin-pair correlation, they are defined as follows [1] 



q k {t) =< <T fc (t) >= J2 a kP(W},t), 



(8) 



7*1 (*) =< <Jk{t)ai{t) >=^a k aiP{{a},t). 



(9) 



According to the definitions and the master equation ([fh, and using the normalized condition (^), time-evolving 

equations of q k (t) and "fki (t) can be derived ( see Appendix 



j t q k {t) = -q k (t) + z2[/2 a kW k {a k -> a k ) J P({a}, t), 

{a} \ & k J 



(10) 



dt 



7h(*)=-27«(*)+X! 



M 



-o\ y^gfeWfc(cr fe -> (T fc ) 



c/ -> cr; 



p(M,t). 



(ii) 



Similarly, time-evolving equation of equal-time multi-spin correlation function can be further derived 
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— < <7i 1 {t)(Ji2{t) ■ ■ -<7 in {t) >= —Tl < Oi 1 Oi^...O in > 

P({a},t). (12) 

Eqs. ([!]), and (|Ic|)~(|l2]) are the basic formulas of the single-spin transition-type critical dynamics which suits 

one-dimensional arbitrary-spin systems. 

All of these formulas can readily extended to spin systems on square lattice and cubic lattice. The only corrections 
are changing subscript i into ij and ijk respectively. 



EE 

M U=i 




III. EXACT RESULTS 



In this section, three examples of the application are given, including one-, two-, and three-dimensional kinetic 
Gaussian model. 

Now we treat the kinetic Gaussian model. First of all we will introduce Gaussian model, then we will exactly 
solve the evolution of the local magnetization and equal-time spin-pair correlation function, and we will obtain the 
dynamical exponent z. Here, we will only give the solving process of three-dimensional case in detail. 

The Gaussian model, proposed by T.H. Berlin and M.Kac [^0| at first in order to make Ising model more tractable, 
is an continuous-spin model. Comparing with the Ising model, besides having the same Hamiltonian form (three- 
dimensional case) 



JY 



ij,k—l w 



(13) 



where meas summation over near neighbors, there are two extension: First, the spin aijk can take any real value 
between (— oo, +00). Secondly, to prevent from all spins tend to infinity, a probability of finding a given spin between 
a^k and a^k + daijk is assumed to be the Gaussian-type distribution 



f(&ijk)d(Jijk 



exp 



'ijk 



d<Jijk, 



(14) 



where b is a distribution constant independent of temperature. Although being an extension of Ising model, Gaus- 
sian model shows much difference from Ising model in the properties of phase transition. In equilibrium case, on 
translational invariant lattices the Gaussian model was exactly solvable, and later as a starting point to study 
the unsolvable models it was also investigated with mean field theory and momentum-space renormalization-group 
method[21-22]. Recently the Gaussian model on fractal lattices was studied by Li and Yang[23]. However, the critical 
dynamic problem of continuous-spin model have never been investigated so far. 

We now proceed to treat the isotropic kinetic Gaussian model on the cube lattice. The system Hamiltonian and 
the spin distribution probability are (15|) and ( |l4| ) respectively. In this case the spin transition probability can be 
expressed as 



W, 



1 



ijk{<?ijk -> Pijkj 



Q 



ijk 



■ exp 



ijk 



ipi+w,j.k ~t~ &i,j-\-w,k ~t~ &i,j,k-\-w) 



Because the spin take continuous value, the summation for spin value turns into the integration 



E 



f(a)da, 



(15) 



(16) 



then the normalized factor Qijk can be determined as 



Q 



ijk 



/oo 
daijkf(aijk)ex.p 
-00 



k&ijk ^ (&i-{-w ,j.k &i.j-\-w ,k ~b &i,j,k-\-i) 



6XP<1 26 



-1 2 



E 



^i-^w ,j.k H~ &i.j-\-w ,k ~i~ ,k-{-w^) 



(17) 
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and the another useful combination formula can also be obtained 



/oo 
&ijkWijk(<Jijk — » &ijk)f(6'ijk)d<Jijk 
- oo 



i ^ ^ (^"i+«/,j',fc ~l~ @i,j-\-w ,k ~T" ,k-{-w ) • 



(18) 



Substituting ( Jlq ) into the following time-evolving equations of the local magnetization and the equal-time spin-pair 
correlation function 



d 

df. 



Qijk(t) = -Qijkit) + ^2 y^AijkWijk(<Tijk -> °i 3 k) ] P({<r},t), 



(19) 



we get 



rff 



Jijk;i'j'k'(t) — ~^Jikj\i'j'k'(t) + E 



0"ij'Jfe I E (ri'j'k'Wi'j'k'i&i'j'k' &i'j'k' 



+ < J i'j'k' I ^^^ijkWijk{^ijk — > 



3I^y fc W — ~Qijk(t) + T (<li+w,j,k + Qi,j+w,k + 1i,j,k+w) ■ 



(20) 



(21) 



7ijk;i'j'k'(t) — —2jijk;i'j'k'{t) + — E (7iife;i'+tuj',fc'(^) + lijk-i' ,j' +w ,k' {t) 



(It 



'Jijk^i' ,j' ,k'-\-w (^) r Yi-\-w,j,k',i / j / k / (^) ~l~ ^i.j+w ,k;i' j' k' (t) ~r~ r yi,j,k-\-w:i' j' k' (0) ' 



(22) 



In order to solve Eqs.(pl|) and (^2|) in the nearest-neighbor interaction case (iu = ±1), we introduce two generating 
functions [1]: 



and 



Fi(A 1 ,A a ,A 8 ,t)= X! AlA J 2 Af^ fc (t), 

i : j,k— — oo 



F 2 (Ai, • • • , A 6 ,t) — ^ A1AIA3A4 A5 Ag r fijk;i'j'k'(t), 



then the Eqs.([H]) and (p3) turn into the following equations with resect to iq and F 2 , respectively, 



-7r-Fi(Ai, A 2 , A3, t) 
at 



■1+ E^+v 1 ) 



■Fi (Ai, A 2 , A3, i), 



(23) 



(24) 



(25) 



P 2 (Ai, • • • , A 6 , i) 



i=i 



F 2 (Ai, • • • , A6, t). 



Solving Eqs.(B5j) and (|26|), we get 

Fi(Ai, A 2 , A 3 , t) = Fx(Ai, A 2 , A 3 , 0)e-* exp 



■E(*+*r x )* 



(26) 



(27) 



F 2 (Ai, ■ ■ ■ , Aa,t) = F 2 (X 1 , ■ ■ ■ , A 6 ,0) e - 2t exp - £ (A, + Ar 1 ) * 
In terms of a generating function of the Bessel functions of imaginary argument 

oc 

e x(A+A- 1 )/2 = £ A ^ (x)j 
^ — — oo 

we obtain immediately the following exact solutions 

-t f 2k \ r f 2k \ r ( 2fc \ 

qijk{t) = e 2^ ?nmi (0)7i_„ I —t J Ij- m I yi I 7 fe _j I yt I 



(28) 



(29) 



(30) 



1ijk;i' j' k' {p) — e ^ ^ 'Jnml-.n'm'l' (0) 



n.m.L\n ,m 



2k 



2fc \ /2fc \ /2fc \ f2k\ (2k 

, i yt I /fc-i I y* 1 W I y* ) I yt I ( y* 



where g„ m i(0) and J n ml;n'm'l' (0), respectively, correspond to their initial values. 

By using of the asymptotic expansion expression of the first-kind imaginary argument Bessel function 



I v (x) 



^TTX ^ n (2x) n 



l^x ^ (2a:) n ' 



where 



(v,n) = 



T^+v + n) 



nT (i + v - n) ' 

we can obtain the long-time asymptotic behavior of the local magnetization 



n. , m . I — — oo 



(31) 



(32) 



(33) 



T = 



l-6fc/&' 



(34) 



where r is the relaxation time of the system. We know that the critical point of the Gaussian model is k c = J/kpT c = 
b/2d, where d is the spatial dimension, and the correlation length critical exponent is v = 1/2 [22]. So, by means of 
the following dynamical scaling hypotheses 



IT - T r . 



(35) 



the dynamic critical exponent z of the 3D kinetic Gaussian model can be obtained 

z = 2. 



(36) 



(37) 



With the same treatment, we can easily solve one- and two-dimensional kinetic Gaussian model. Ignoring the 
process of solution, we give only the following exact results: 
(1). ID case 



G 



~ 2k 

q k (t) = e~ t 22 5m(0)I fc - m (yt), (38) 



>0 



2A; 2k 

lu{t)=e- 2t V 7„ m (0)I fe _„(— f)/j_ m (— i), (39) 



y , I i bi 1 1 \ " / ~ rv — i b \ j ) I — lib \ j 

b b 

n.m— — oo 



T = T^~2kjb' (40) 
z = 2. (41) 



(2). 2D case 



°° /2fc \ / '2k \ 

qnm(t) = e- t Qk,i(0)I n -k \^-tjI m -i i^—tj , (42) 

7mn;m'n'(tJ = e 2^ 7ij;i'i'(°J J m-i I "7"* I Ai-j I "7"* I hn'-i' I "7"* J i n '-j' I "7"* ) (43) 

,j'=— oo 

T = T^> (44) 

z = 2. (45) 
IV. CONCLUSION 

In this paper, we have suggested a normalized transition probability to solve the time evolution equations of the 
local magnetization and spin-pair correlation function. Our treatment can in part remove the arbitrariness in Glauber 
dynamical theory, and makes exactly solve the time evolution equation possible. 

Based on our treatment, we have re-investigated one-dimensional kinetic Ising model, resulting in same dynamical 
critical exponent z as Glauber[l]. We have also exactly solved the kinetic Gaussian model, and given the detail 
of solving process of three-dimensional case. We have in surprise found that the dynamical critical exponents have 
same value independent of spatial dimension, which shows the dynamical behavior has super-universality in con- 
trast with static behavior. In fact, in equilibrium phase transition the critical exponents are strongly dependent on 
dimensionality. 
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APPENDIX A: PROOFS OF EQS. (10)-(12) 

According to the definition (B) and using the master equation (11]), we have 
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M 

{cr} * cr; 

= £°'fc £ / ] [-Wi(a-j -> o-j)P({o-}, t) + Wi(ai — > g-f )-P({q>v»}, dj, *)] 

{cr} i{i=£k) (Ji 

+ [ _<T & U, &( <J fe 0-fe)P({cr},i) +(TfeWfc((7fc -> CTfc)P({0"j^k},0-fc,t)] (Al) 

{cr} rTfc 

Looking at the first term (i 7^ fc) in the latter of the last equality sign of flAl| ), 

(i 7^ fc) term = £o-fc £ ^[-^(gj -> f7i)P({cr},t) + u>i(<7i -> CT i )F({cr i ^ i }, 5*, t)] 

{cr} i(i^k) ai 

= £ CTfc £ [- £ Wi(o-i -> a J )P({cr},<) + £ Wi(<7; a^P^a^},^,*)] 
{<Jj^i\ i(i^k) <yi,<ii Oi.Hi 

(A2) 

it is easy to see that this term equals to zero, as long as &i exchange with Gi before doing sum for <7j and Oi. So the 
surplus term of (Al) is only the last term (i = k): 

= [~ a kW k {ok —> a k )P({cr},t) +(J k w k (& k — > (Tk)P({(Tj3&k},<rk,t)] 

{cr} cr k 

= - yVfc 1 y^^fc(gfc -> 1 + ^kw^k— ^^({0^},^^) 

{cr} \ cr fc / (71 ■■■<J k cr k ■■■cr N 

= -y^gfcP({cr},t) + ^ OkWk(<Jk -> o'fc)P({o'i#fc},o'fc,t) 



{cr} a\---a k o k ---ON 

-Qk( 



{cr} V <Tfc / 



{cr} \ cr fc 

in which, the normalized condition and the technic of exchange of <7i for cr^ were used. Hitherto, Eq.(lO) have been 
proven exactly. As for the proof of the Eqs.(ll) and (12), it is easily accessible via the same method, so without the 
necessity for further proof. 



[1] R. J. Glauber, J. Math. Phys. 4, 294 (1963). 
[2] K. Kawasaki, Phys. Rev. 145, 224 (1965). 

[3] G. Forgacs, S. T. Chui and H. L. Frisch, Phys. Rev. B 22, 415 (1980). 

[4] E. J. Lage, J. Phys. A: Math. Gen. 18, 2289 (1985). 

[5] E. J. Lage, J. Phys. A: Math. Gen. 18, 2411 (1985). 

[6] E. J. Lage, Phys. Lett. A 127, 9 (1988). 

[7] Y. Achiam, Phys. Rev. B 31, 4732 (1985). 

[8] Y. Achiam, Phys. Rev. B 32, 1796 (1985). 

[9] Y. Achiam, Phys. Rev. B 33, 7762 (1986). 
[10] Jian Zhou and Z. R. Yang, Phys. Rev. B 39, 9423 (1989). 
[11] M. D. Lacasse, J. Vinals and M.G rant, Phys. Rev. B 47, 5646 (1993). 
[12] B. C. S. Grandi and W. Figueiredo, Phys. Rev. E 54, 4722 (1996). 



8 



[13] J. Rogiers and J. 0. Indeleu, Phys. Rev. B 41, 6998 (1990). 
[14] J. Wang, Phys. Rev. B 47, 896 (1993). 

[15] Y. Achiam and J. M. Kosterlitz, Phys. Rev. Lett. 41, 128 (1978). 

[16] G. F. Mazenko and O. T. Vails, Phys. Rev. B 24, 1419 (1981). 

[17] P. C. Hoohenberg and B. I. Halperin, Rev. Mod. Phys. 49, 435 (1977). 

[18] Z. R. Yang, Phys. Rev. B 46, 11578 (1992). 

[19] Peiqing Tong, Phys. Rev. E 56, 1371 (1997). 

[20] T. H. Berlin and M. Kac, Phys. Rev. 86, 821 (1952). 

[21] J. J. Binney, N. J. Dowrick, A. J. Fisher and M. E. J. newman, The Theory of Critical Phenomena, Clarendon Press, 
Oxford 1992. 

[22] L. E. Riechl, A Modern course in Statistical Physics, University of Texas Press, Austin, 1980. 
[23] Song Li and Z. R. Yang, Phys. Rev. E 55, 6656 (1997). 



9 



